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We consider a model of spin-gapped chains weakly coupled by Josephson and Coulomb inter- 
' actions. Combining such non-perturbative methods as bosonization and Bethe ansatz to treat the 

, intra-chain interactions with the Random Phase Approximation for the inter-chain couplings and 

the first corrections to this, we investigate the phase diagram of this model. The phase diagram 
shows both charge density wave ordering and superconductivity. These phases are separated by line 
of critical points which exhibits an approximate an SU(2) symmetry. We consider the effects of a 
magnetic field on the system. We apply the theory to the material Sr2Cai2Cu24 04i and suggest 
further experiments. 
PACS numbers: 71.10.Hf, 74.20.Mn 



I. INTRODUCTION 



^ . ..... 

Quasi-one-dimensional (ID) models are often used to test various theoretical ideas in the area of strongly correlate] 
. electron systems for the simple reason that most known non-perturbative approaches work only in one dimension.GI'1 
The route often taken is to use a non-perturbative solution of a strictly one-dimensional model and then use mean field 
r , or the Random Phase Approximation (RPA) to take into account the inter-chain interactions. Through techniques 
' ^ ' such as bosonization and Bethe ansatz, many results are known about such one-dimensional systems as spin chains 
and Tomonaga-Luttinger liquids which form the skeleton of these quasi-one-dimensional models. Linking thesCr-using 
the RPA formalism hafi-yielded many successful experimental predictions, for example for linear conductors and 
for magnetic systemsJflO A very early use of this technique is Efetov and LarkiriEla who estimated the transition 
temperatures in the same model we use. 
' As it is well known, the RPA formally represents the leading term in a perturbation expansion in where z±_ 

, is the number of nearest-neighbour chains in the lattice. For real experimental systems this number is not usually 
\^ ' large so it is important to know about higher order contributions in \/z^. The recent results for the quasi-one- 
^NJ , dimensional tJeisenberg magnets indicate that the worst these corrections can do is about a 25% shift in the transition 
T-H I temperatureBM The relative smallness of these corrections demonstrates the validity of the RPA approximation when 
considering real systems; in our case it turns out that the corrections are even smaller. 

In this paper we follow the same road and discuss a simple model of a non-BCS superconductor. In the model we 
consider the formation of superconducting pairs on one-dimcnsional chains is triggered by formation of a spin gap. 
The three-dimensional coherence is established through the inter-chain Josephson coupling. Since the latter coupling 
competes with the Coulomb interaction, which can destroy the superconductivity and establish Charge Density Wave 
(CDW) ordering. As we shall show, these two phases are separated by a critical line with increased symmetry. Near 
this line, we take into account the interplay between these two interactions considering corrections to RPA. 
P3 The model we use has been considered in some detail recentlyliJ in the context of higltTc superconductivity. It was 

O ■ assumed that the one-dimensional behaviour came abO|U± from the formation of stripes.tHI Since in the stripe picture, 
^ ' fluctuations of the stripes dephase the CDW couplingJi3 only the SC inter-chain interaction was considered. In our 
^ . paper we retain the Coulomb interaction and therefore expect it to be relevant to materials that are structurally 
' quasi-one-dimensional such as the Bechgaard salts or some cuprate materials such as the family Sri4_xCaxCu2404i. 
rS , In section || we introduce the model we will be dealing with. In section III we show that this model has an SU(2) 
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' symmetric quantum critical line. In section tV we calculate the transition temperature for the model within the RPA 
approximation. Treating the inter-chain coupling in the mean field approximation we obtain an effective sine-Gordon 
model for each chain. Using the exact results for this model we calculate the zero temperature spectral gap M and 
derive the expression for the ratio Tc/M. Here, we also consider the properties of our system in a magnetic field. 
In section ^ we look at the first corrections to RPA which gives us an improved phase diagram of the model. In 
section ^ we show that the same general behaviour also occurs in two dimensions, although the transition here 
is Kosterlitz-Thouless rather than the symmetry breaking found in higher dimensions. Finally, in section VII, we 
show that the quasi-lD compound Sr2Cai2Cu2404i is a beautiful example of our model and we discuss the measured 
properties of it in relation to our theory. We also make some quantitative predictions about this material which could 
be confirmed by further experiments. 
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II. THE MODEL 



Let us consider a system of conducting one-dimensional units (we will call them chains, though in reality they may 
be, for instance, ladders) weakly coupled to each other. As often happens in one dimension, the spin and charge 
degrees of freedom decouple at low energies. We assume that the spin sector aquires a gap, and that the filling of each 
individual chain is incommensurate with the lattice so that the low energy behaviour of the charge sector is decribed 
universally by the Gaussian model. The Hamiltonian density is therefore 

'^chain — ^charge ^" '^spin; (1) 
^charge = ^[K,{d,ef + R-^d^^f] (2) 

where [Q{x),^{y)] = i9{x — y). We don't write down the form of the Hamiltonian in the spin sector: our only 
requirement is that it has a gap A^. [] 

The spin gap blocks single-particle tunneling processes between the chains. Then the multi-particle processes 
generate pair hopping. In what follows we shall assume that the inter-chain tunneling matrix element is much smaller 
than the spin gap. In this case one can take into account only two-particle virtual processes giving rise to Josephson 
coupling between the chains. They lead to the following Hamiltonian 

T-isc = \ ^ : cos[\/27r(e„ - Q,n) - 2eHbnmx/c] : (4) 

where the dots signify that operators are normal ordered with respect to the state with spin gap and therefore the 
ultraviolet cut-off for the correlation functions of bosonic exponents is Aj. The fields without index are assumed to be 
from the charge sector, as will be the case from here on. We have also introduced external magnetic field H directed 
perpendicular to the chains; 6„m is the projection of the inter-chain lattice vector on the direction perpendicular both 
to the chains and the magnetic field. 

An analysis of dimensionalities as shown in Appendix ^ yields 

A,\^/^^-^ 



^-^Uj a: 

where t is the single particle hopping and A is related to the original bandwidth. 
Interaction has scaling dimension 

dsc = l/(2i^c) (6) 

and therefore is relevant even for repulsive interactions in the charge sector provided they are not too strong [Kc pa. 
1/2). This is a well known effect of the spin gap; it generates preformed pairs making it easy for them to condense.Efl 
There is also a Coulomb interaction between the two chains In the spin gap regime, there is only one term in here 
that remains relevant: it is the coupling of 2A:i;--components of the charge density which gives the effective Hamiltonian 
density (Appendix ^) 

Hcdw = ^Kff J2 ■ cos[V2^($„ - $™)] : (7) 
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where 



^In a specific case of single chains a realistic description of the spin sector is given by the SU(2) Thirring model Hamiltonian :EJ 

H.pin = ^(: rr : + : :) - : J" J" : (3) 

where Vs is the spin velocity and J'^,J°' are chiral SU(2) currents satisfying the level k = 1 SU(2) Kac-Moody algebra. The 
spin gap is generated when g > such that the current-current interaction in the spin sector is marginally relevant. In the 
case of ladders a description of the spin sector is more complicated; this however does not affect the charge Hamiltonian and 
therefore will not concern us here. 
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-y-f) ^0. (8) 

The corresponding scaling dimension is 

rfcdw = KJ2. (9) 
The effective action for coupled chains is therefore 

£eff = ^[^^,^nf + \Y^ {Vnm ■ cos[V2^($„ - : + J„™ : cos[V2^(e„ -Sm- 2eHhnmx/c)] :} (10) 



2K^ ' 2 

and has As as the ultraviolet cut-off. We will be considering nearest-chain interactions only, i.e. Vnm — V, Jnm — J 
for neighbouring chains and zero otherwise. In what follows we will be most interested in the case Kc ~ 1 when both 
interactions are important. 

III. AN EFFECTIVE THEORY OF THE CRITICAL POINT 

For a general value of Kc the symmetry of the model is U(l)xU(l) which corresponds to independent global 
shifts of <& and 8. When Kc = 1 and ±J the symmetry increases and becomes SU(2). To see this we use 

the non-Abelian bosonization descriptionBt^ At Kc = 1 the exponents exp[±i'\/27r$], exp[±i-\/27r0] have conformal 
dimensions (1/4,1/4) and can be understood as matrix elements of the tensor field gab from the S=l/2 representation 
- the first primary field, of the level k = 1 Wess-Zumino-Novikov-Witten model (for a discussion of this model, see e.g. 
Itzykson and DrouffeEEl): 

„_/ exp[iV27r$] exp[iV27re] \ , , 



exp[-iv27r0] exp[-iv27r$] 
The Gaussian part of the action becomes the sum of the WZNW actions from individual chains: 

n n 

and the interaction term in ( p^ can be written as 

Unt = E ~ E [5™]^""^ + (n ^ rn)\ + JTr(5„g+ + g^g^)}. (13) 

n^m a— 1,2 

This description is convenient since it contains only mutually local fields and therefore can be considered as the 
Ginzburg-Landau theory. 

In three spatial dimensions the system undergoes a phase transition into the ordered state where the matrix g 
acquires an average value throughout the system. In the long wave limit one can replace the last term in (|l3| ) by 

{dyg){dyg+) (14) 

and omitting the time dependence of the fields we obtain the following Ginzburg-Landau free energy: 

P = ^'^ j da;dVTr[^(a,<?+9,5) + Jb\V ^g+V ^g)] + Famsotropy (15) 

where h is the lattice constant in the transverse direction and 



Fanisotropy 



{V-J)b-'^ I dxdV ^ g^'"'\g+]'^'"'\ (16) 



a=l,2 



We can now re-parameterize the theory. The order parameter is the SU(2) matrix g. Its relation to the GDW and 
SC phases 8 and $ are: 

g = exp[icr3($ + 8)/4] exp[icria/2] eyij>[ia^{^ - 8)/4]. (17) 
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The Ginzburg-Landau free energy density is 

T ^ ip[cos2(a/2)(Ve)2 + sin2(a/2)(V$)2] + ^^(Va)^ + {V - J) cos, a. (18) 

This is interpreted as foUows: when y — J is positive, a is pinned at tt so that the coefheient in front of (V<i>)^ is 
non-zero and hence $, the CDW order parameter, is constant throughout the material. When V — J is negative, a 
is pinned at and hence it is 0, the superconducting order parameter that acquires an expectation value. When 
y — J = we are at the critical point where the free energy of the superconducting and insulating phases becomes 
equal. The effects of this V — J mode will be considered throughout the rest of the paper. 



IV. PHASE DIAGRAM IN MAGNETIC FIELD AND CRITICAL TEMPERATURE 

For two chains the problem was solved by Shelton et aZ0 There are two modes; one symmetric in the two chains 
and the other antisymmetric. In the presence of the inter-chain interactions, the symmetric mode remains gapless 
and the antisymmetric sector splits into two Majorana fermions with gaps (V -f J) and (V - J). 

For an infinite number of chains, we expect to see a similar sort of behaviour. The gapless symmetric mode in the 
the case of two chains will in some sense be the Goldstone mode in our infinite system and we expect to see a range 
of other modes with gaps ranging from V — J \oV ^ J . We will see that within the basic RPA approximation we 
cannot reproduce this behaviour: the properties will depend on the stronger of V and J but not both. However when 
we go beyond the first order term we can start probing the interplay between these two competing interactions. 

To begin with, we estimate the critical temperature using RPA. Within this approximation the pairing and the 
CDW susceptibilities are given by 



Xsc 



1 - Jz_Lxi? 



^(0) 

Xcdiu = TJyp (19) 

where zj^ is the number of nearest neighbour chains. These are shown diagrammatically in figure ^ (a) and (b). 

When Kc = 1 the bare susceptibilities are equal to each other and therefore the instability occurs in that channel 
where the interaction is stronger. This is shown explicitly in Appendix^. If Kc ^ 1, the instability still occurs in the 
stronger channel, although this now depends not only on the values of V and J but also on and , the crossover 
point being 

An important modification occurs in magnetic field which affects the inter-chain interaction in the superconducting 
channel (Q). In this case the susceptibilities corresponding to the lattice directions 1 should be taken at wave vector 
2e(H[x X l])/c, where x is the unit vector along the chains. Therefore the RPA criterion for the transition is replaced 

by 

l = E^'^-{9 = 2e(H[xxl])/c} (21) 
/ 

For definiteness let us assume that the instability occurs in the superconducting channel which is the most likely 
case for Kc > 1. Note that the duality property of the effective Lagrangian (|l^) under K ^ 1/K, V ^ J, 6 ^ $ 
means that all of the results in this and the next section are identical for the CDW channel. 

In a Tomonaga-Luttipger liquid with the ultraviolet cut-off Ag the static susceptibility for the operator with scaling 
dimension d is given byll3 



X(°)(<z) = -^sin7rd('^') T\l-d) 



^^^'^ Tid/2 + ivq/ATTT) ^ 



where v is the velocity in the charge sector 



T{l-d/2 + ivq/ATTT) 



(22) 
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A. Zero magnetic field; tiie critical temperature and tlie vortex energy 

Substituting ( p^ with q — into Eq.(p^) we obtain 

^ A, /2jzi . r2(d/2)r2(i-d) 

J^c — t: — — : — sinvra — —77 -7— 

27r V As r2(l~d/2) 



(23) 



The scahng properties of this equation were calculated for the first time in reference ^ 

Below the transition temperature the long-wavelength fluctuations of superconducting order parameter are three- 
dimensional. The amplitude fluctuations are, however, mostly one-dimensional and their spectral weight is concen- 
trated above certain energy which plays a role of a pseudo-gap. The zero temperature value of the pseudo-gap can be 
found from the mean-field theory combined with the exact results for the sine-Gordon model. In this approach one 
approximates the inter-chain interaction 



J ^ COS^(0„ - (j),n) 



<nm> 



( (3(j) = V2^e and /J^ = 2nK-^) by 



where 



2/x cos / 



2/^ = J z± As{cos (3(j)). 



This expectation value is known exactljJil: 

(i + 07rr(i™d/2) /r(i + |)r(i-|; 



(d-2) 



(cos /3</)) 



16sin< r(d/2) 



2sin- 



atJ 



(24) 

(25) 
(26) 

(27) 



where M is the soliton mass in the SG model, and is related to /i by 

d-2 



^ r(d/2) / 2r(e/2) 
^ 7rr(i-d/2) 1^0Fr(i + |) 



4^V"a?. 



(28) 



In all these equations, d = 0^ /Att is the scaling dimension of the field e^^"^ , and ^ — f/(2 — li). These mean- field 
relations are solved to give 



M = A, 



JZ± 1 TT^ 

— -, tan — 

As 2(fi- 2) 2 



7rr(i-d/2) /r(i + |)V^' 



{d-2) 



V{d/2) 



2r(e/2) 



(29) 



The ratio Tc/M which is often considered in the theory of superconductivity is plotted as a function of d in figure 
|l|. It's numerical value in certain limits is: 



§(..0).f .0.177, 

T^,^_ ,/2) - 3 V3^(r(2/3)r(5/6))3 
A/^ /^"16 r(3/4)8 



(30) 
(31) 



In the limit d — s- 1 which corresponds to weak coupling, our expressions for and M diverge in this approximation. 
However their ratio can still be evaluated. Writing x ~ \ — d and expanding all the gamma functions as Taylor series 
in X gives us the BCS value 



§(d ^ 1) = hm [1 + (ln2 + 7)0;]* = « 0.567 



(32) 



where 7 ~ 0.57722 is Euler's constant. 
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Notice that in comparing to experiments, one has to remember that M is not the single particle gap. Single particle 
spectroscopies such as basic tunneling would see a gap much closer to Ag, the spin gap. To probe M, one would 
have to look at experiments involving pairs of electrons, such as Andreev tunneling. In the context of sine-Gordon 
model, M is the soliton mass. Solitons correspond to spatial changes in the superconducting phase O and hence to 
vortices. Therefore M is the minimal energy necessary to create a vortex. It should also be noticed that at d < 1 
the sine-Gordon model has not only solitons, but bound states which, being neutral, should be interpreted asAwrtex- 
antivortex pairs. At c? < 1/2 the energy of the first bound state is smaller than the soliton. See Carlson et aZO for a 
nice discussion of the implications of having these two energy scales. 




0.2 0.4 d 0.6 0.8 



FIG. 1. A graph of Tc/M against d. The value d = 1 corresponds to the BCS limit, decreasing d corresponds to increasing 
repulsion. 



B. Phase diagram in a magnetic field 



To keep the calculations as simple as possible, let us consider the simplest possible situation when a given chain has 
four nearest neighbours with Josephson couplings Jz and Jy and the magnetic field lays in the yz plane. Combining 
Eq.(pT[) and Eq.(p^) we obtain the equation for the critical temperature: 



C 



{2-2d) 



= Jz 



C 



T{d/2 + tabzHy/T,) 



r(l -d/2 

= (Jz + Jy 



iabzHy/Tc 

r(d/2) 



r(i - d/2) 



Jy 



a 



T{d/2 + iabyHz/Tc) 



r(l-rf/2 
: ew/27rc 



iabyHz/Tc) 



(33) 



The solution of this equation describes several interesting effects. 

• A possibility of a re-entrance behaviour. 
Let us consider the case when in-plane interactions are isotropic: 



Jz 



Jy:^z 



by and the magnetic field is 



directed at 45° angle — Hy = H . This gives it the maximal power to suppress T^. A numerical solution of 
Eq.(^3|) is plotted in figure || (a) for various values of the scaling dimension d. We see that there is a range of 
magnetic fields for which the superconductivity exists in an intermediate range of temperatures. To study the 
stability of these solutions one needs to have a good description of the ordered state in magnetic field, which we 
hope to obtain in the future. 

At Tc ^ Eq.(p3|) can be solved analytically which allows us to extract the value of critical field at Tc = 0: 



2TrcTc{0) /r(l -d/2)\ 
~e hT \ r{d/2) J 



i/{i-d) 



(34) 
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This is plotted in figure || (b) along with the numerical solution for 
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FIG. 2. (a) The critical temperature as a function of magnetic field for various values of d. (b) The critical magnetic field as 
a function of d. The magnetic field is measured in the units of 2ehbv / c. 



• Anisotropy of the phase diagram. 

Another prediction following from Eq.(p3|) is an anisotropy of the phase diagram. This can be illustrated by an 
analytical solution for Tc ^ case. Setting Tc — > Eq.(p3|) in we find 



J. 



J II 



C 



[aHyh.yi^-d) {aH,by)^(^-d) [r,(0)]2(i-'^) 



(35) 



This is plotted in figure ||. We must be careful to remember however that this is a first order mean field 
calculation, and further corrections will give a critical flux in all directions, even when the field is pointing 
directly along one of the crystal axis. 

• SC-CDW transition. 

The validity of the above calculations is limited by the range of temperatures where the system is stable against 
CDW transition. Therefore, strictly speaking, before the Tc quantum critical point is reached the system 
will undergo a transition into a CDW state. 



FIG. 3. Angular dependence of the critical magnetic field, 
other values of d. 



This is plotted for d = 1/2. The graph is qualitatively similar for 
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V. CORRECTIONS TO RPA 



The analysis of the previous sections was based on RPA. Since in reahstic situations the number of nearest neighbours 
is never large, it is important to check how robust RPA is. We will calculate corrections to RPA in the simplest case 
case of zero magnetic field. We shall also restrict ourselves to Kc = 1 (d = 1/2 for both interactions). 

The basic RPA calculation involves only the stronger of the two interactions - for clarity let us again take this to 
be J. However as we mentioned before we would expect the presence of the other competing interaction of the same 
scaling dimension to also play a role. In particular we expect there to be a mode with a gap of J — V, seen in ( p^ ) 
and in the two chain model. This will be very important around the point y = J as it will become massless thereby 
increasing fluctuations and decreasing the transition temperature. This can be investigated by looking at the first 
correction to the RPA formula - figure 0(c). 



(a) 




+ 



(b) 




(c) 




+ 



FIG. 4. (a) The basic RPA diagram, (b) The Dyson series for RPA, (c) The first correction term. In these diagrams, the 
dashed lines represent the ID chains, the dots indicate vertex operators of or and the wiggly lines are the inter-chain 
interactions, and each diagram is an irreducible correlator. 



In terms of the fields </> and 9, this diagram can be expressed as 

5x - V^z^ 



^iV2^<f,(a) ^iV2^9(l) ^-iV2^9{2) ^-iV2^4>{b)\ _ ,iV2^4,{a) ^-iV2^4>{b)\ , 1^2^9(1) ^-iy^ 8(2) \ 



^^g-iV2¥e(l)gjV2¥9(2)^ 



The revised RPA equation for the transition temperature is 



X (e 



-iV27r0(l)giV27r0(2)\ 



(36) 



(37) 



where the coefficients are given by 

dr / dx- — — — 

/o ' 



^0 



1 /"^ /"OO 

— / dr I dx — , M „ 

TT .In J_oo |smh(a; + ir)| 27r 



=.^52(1/4,1/2) 



1 







dTidT2dTi, / dxxdx2dxb 



sinh(a;6 + iTf,)| | sinh(.Ti2 + iTi2)p 
sinh(a;i + iri)| | sinh(xf,2 + in2)\ 



2 /"TT POO 

AY = — dTidT2dTb / dxidxidx, 

Jo J-oo 



I sinh(a;2 + ir2)|| sinh(a::f,i + ini)\ 
1 1 



sinh(a;6 + in)] \ sinh(xi2 + jTi2)p 

1/2 

/ sinh(a:i + iri) sinh(a;2 — «T2) sinh(xb2 + iT62) sinh(a;bi — iTfti)^ 
\sinh(a;i - iTi) sinh(2;2 + ^'''2) sinh(xb2 - ir^) sinh(a;bi + iTbi) 



(38) 



(39) 



(40) 
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with Xi2 = X2 — Xi and so on. _ 

The integrals are evaluated numerically be Monte-Carlo techniqueSjCB with values calculated over finite volumes 
then scaled to infinity. The results are 



A'^ = 4.377, 

A( = 34.81 ±0.02, 

A]^ = -33.01 ±0.02. 



Hence the correction to the transition temperature is 



1 + — 



0.42 - 0.40 — 



(41) 



(42) 



This expression is valid for J > V . liV>J, the expression is exactly the same, but with V ^ J . This is plotted in 
figure H and gives a dip near the critical point as expected. 

It is interesting to note that in the absence of the second interaction term, ie = 0, these correction raise the 
transition temperature above the RPA Yft^*^' This differs from models of coupled spin chains where RPA tends to 
overestimate the transition temperaturcEla 




0.1 



0.05 



0.75 - 




0.25 



0.05 
V 



CDW 



0.1 



FIG. 5. (a) A plot of Tc against V and J, (b) A cross section of Tc against V along the line V ~\- J = 0.1. In these plots, we 
have taken zx = 2 to allow these corrections to be clearly seen, although for this approach to be valid, we require z±_ > 3. 



VI. A WORD ABOUT TWO DIMENSIONS 



In two dimensions the RPA approach in the previous two sections must break down completely, as spontaneous 
symmetry breaking is forbidden by the Mermin- Wagner theorem. We can see how this comes about by looking at 
figure ^(c) . The correction we looked at involved only bare couplings to the bare correlation function. The process 
of making these lines 'thick' involves much numerical complication and gives rise to only small corrections in three or 
higher dimensions.^ However in two dimensions, these corrections have infra-red divergences and drive the transition 
temperature back down to 0. nn 

Nevertheless we still get a transition in two dimensions: it is of the Kosterlitz-ThoulessEllES type. Let's look closer 
at Coulomb coupling in two dimensions. The Lagrangian for the coupled chains can be written 

^ = E {^(^M</'0' - Jcos[/3(<^, - (/),+i)]| . (43) 
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By making the approximation 



cos V — ~2 ' 



(44) 



which comes from the diagrammatic expansion, we can write this as 



(45) 



with the self-consistent relation 

J = J/32 (cos /3(,^,- 0^+1 )) 



1 — COS q±^ 



Lol + ql + AJsm\q^/2) 



(46) 



At T = this relation becomes 




J = J/J^exp -^In 



J/32 



'2J , 



(47) 



where d = 0^ /Att as before. As T increases, the self-consistent value of J will decrease, but for an estimate of the 
behaviour of the transition temperature this relation will suffice. 

The Kosterlitz-Thouless transition temperatureEll^l Tkt ~ \/j hence we have 



Tkt - A, 



J_ 

a: 



(48) 



giving the same order of magnitude as the ordering temperature in higher dimensions ( p3| ) . 

Hence in two dimensions, although the nature of the transition is different, the energy scales involved are the same 
as in higher dimensions. The only major difference occurs when approaching the SU(2) critical point where the 
presence of a non-Abelian symmetry in two dimensions means that the transition temperature will drop to zero at 
this point. The qualitative phase diagram in two dimensions is shown in figure ^. 




V/J 



FIG. 6. The modified phase diagram for our model in two dimensions 



VII. AN EXAMPLE EXPERIMENTAL SYSTEMS 



The class of materials Sri2-xCaxCu2404i are built up from alternating layers of weakly coupiejd Cu02 chains and 
CU2O3 two- leg ladders. The material shows a spin gap in both of these one-dimensional units,E2l making it a prime 
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candidate for application of our model. Our theory is still valid even if the superconductivity originates from the 
ladders. r~ir~| ri 

For X > 11.5, these materials show superconductivity under pressure, cilEa and NMRE3 studies also indicate possibkj 
charge ordering at low temperature and ambient pressure. Recent measurements of the electrodynamic responsecJ 
have confirmed the presence of CDW in this class of compounds. One of the most interesting measurements however 
is the DC resistivity. For Sr2.5Caii.5Cu2404i these measurementsE^ show a number of features: 

• Below about 4 GPa pressure, the temperature dependence of the resistivity perpendicular and parallel to the 
ladders is different. This indicates that different mechanisms are governing the transport in these two directions, 
consistent with the spin-gap concept. Above 4 GPa the temperature dependence of the resistivity anisotropy 
becomes weak, which indicates that single particle hopping between ladders is now possible, i.e. the spin gap 
has vanished and we have a crossover to a conventional two-dimensional metallic bebaviour. This is consistent 
with the pressure dependence of the spin gap observed in recent NMR experiments.EZI 

• At sufRciently high temperatures, coherent inter-ladder charge dynamics is also seen. The temperature where 
this occurs is consistent with the NMR determinations of the spin gap, so we may conclude that the transport 
properties of this material are indeed governed by weakly interacting one-dimensional spin-gapped units. 

In figure |^ a qualitative phase diagram of this material is shown.0 This is explained in terms of our model. If we 
take if s R:! 1 we have 



JcS ^ f^/^s, 

V,s ^ Fo(A,/A). 



(49) 



The increase of spin gap leads to decrease in the effective inter-ladder Josephson coupling. Hence eventually the 
inter-ladder Coulomb interaction takes over and the superconductivity disappears. In quasi-two-dimensional the SC 
and CDW regions of the phase diagram are separated by the quantum critical point, as described in Section VI. 

It would be interesting for this material to measure the charge gap in the superconducting region. This may be 
achieved via optical conductivity measurements. For the Luttinger liquid parameter Kc ~ 1, our model then predicts 
the ratio Tc/Ac to be the non-BCS value of order of 0.4. 

Also in this material, Tc is very small in comparison to the Fermi energy v/a, so the magnetic field effects on the 
superconducting state should be strong. This would be another interesting experiment to perform. 
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FIG. 7. Qualitative phase diagram for the spin gap and superconducting transition temperature against pressure in 
Sr2Cai2Cu2404i 



VIII. CONCLUSION 

We have discussed a model with the following hierarchy of energy scales: 
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1. The highest energy scale is the spin gap Ag. Below Ag the system is described by competing CDW and SC 

fluctuations. 

2. There is a transition temperature at which either (cos v^SttB) or (cos ^/2tt<^) are formed. According to the mean 
field calculation, these order parameters cannot be formed simultaneously. Thus we are either in CDW or SC 
phase, but the temperature of their formation goes smoothly through the point V = J. 

3. There is a third energy scale associated with the gap for another mode which becomes soft at the critical point. 
This mode is not seen in the first order RPA calculations, but it's effects can be noted by looking at the first 
correction to RPA. 

Within the RPA approximation we calculated the transition temperature for general Kc- We calculated the ratio 
Tc/M where M is the zero temperature gap in the charge sector. We saw that this decreases below the BCS value 
as the coupling strength is increased. We also looked at the properties of our model in a magnetic field, noting in 
particular the extreme anisotropy of the phase diagram. 

We then went on to calculate the first corrections to Tc in the vicinity of the critical point which is decreased 
because of the interplay between the two interactions. We also showed that in two dimensions where RPA breaks 
down completely, we get a transition of the Kosterlitz-Thouless type which has the same energy scales as the ordering 
transition in higher dimensions. We also showed that the compound Sr2Cai2Cu2404i is likely to be described by 
our model and on this basis made further predictions about it's properties and suggested that optical conductivity 
experiments should be done on such a material. 

S.T.C. acknowledges many useful discussions with Ralph Werner. A.M.T. is grateful to Marc Bocquet for giving 
him his preprint before publication. The work is supported by the US DOE under contract number DE-AC02-98 CH 
10886 and EPSRC grant number 99307266. 



APPENDIX A: ESTIMATE OF THE EFFECTIVE COUPLINGS 



The easiest case is that of the Coulomb coupling. In the bare system, we have a term 

Wcoulomb = — Pn{x)pm{x) (Al) 

with p{x) the charge density on each chain, and Vq is the strength of the inter-chain Coulomb coupling. When we 
open a spin gap, two things happen to this expression. Firstly, anything involving the spin field is replaced by it's 
average value (cos(\/27r$s)) ^ (As/a(7^)^°/^. Secondly, the cut-off in the normal ordering of the charge sector is 
changed from ag ^ to Ag. This gives an extra factor of (As/ag ^)^'^/^ for each operator. Overall, we generate an 
effective interaction 

Wcdw = \^y,- cos[V2^{^n - : (A2) 
^ As , 

where 

Kff ~ y-^j Vo. (A3) 

In this paper, we will be keeping Kgi^l. 

In the case of the effective Josephson coupling, we start from a single particle hopping term in our bare Hamiltonian 
density 

Whopping = 7, {RnRm + L\iLm) ■ (A4) 

zaQ ^ 

After opening the spin-gap, the effective Hamiltonian density only involves pair hopping: 

Wsc = TT^J^s • cos[V2^(e„ - e„)] : . (A5) 
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These are virtual processes involving an intermediate energy A^, hence the J^s will have a factor t^/Ag. We must 
also remember to change the cut-off in the normal ordering, so the overall expression is 



On 



J.^-\^] ^ (A6) 



APPENDIX B: MEAN FIELD SOLUTION FOR MANY CHAINS 



In the mean field approximation, the interaction term is 

Ant = C0S[V2^(4>„ - + JcOs[V2^(en - O^)]} 

K, ^;_lF(cos[\/2^$]) cos[\/2^$„] + J(sin[\/2^©]) sin[\/2^©„]. (Bl) 

This can be written as 

£i„t = V'l2TKB2Tr[(cos7/ + icr^ sin 7)5 + c.c], 

A= F^_l(cos[\/2^$]), B = j2_L(sin[\/2^e]), tan7 = — (B2) 

The constant matrix can be removed by the redefinition of g. After that it becomes evident that the free energy 
depends only on = + B"^. The mean field equations are 

OF AdF 

^---^'^dB-'-^'^RdR ^^^^ 
Prom this it is clear that the only case where both A and B are simultaneously non-zero is V = J. 
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